Ab Initio Calculations of Even Oxygen Isotopes with Chiral Two- Plus Three-Nucleon Interactions 
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We fonnulate the In-Medium Similarity Renormalization Group (IM-SRG) for open-shell nuclei using a 
multi-reference formalism based on a generalized Wick theorem introduced in quantum chemistry. The result- 
ing multi-reference IM-SRG (MR-IM-SRG) is used to perform the first ab initio study of even oxygen isotopes 
with chiral NN and 3N Hamiltonians, from the proton to the neutron drip lines. We obtain an excellent reproduc- 
tion of experimental ground-state energies with quantified uncertainties, which is validated by results from the 
Importance-Truncated No-Core Shell Model and the Coupled Cluster method. The agreement between concep- 
tually different many-body approaches and experiment highlights the predictive power of current chiral two- and 
three-nucleon interactions, and establishes the MR-IM-SRG as a promising new tool for ab initio calculations 
of medium-mass nuclei far from shell closures. 

PACS numbers: 13.75.Cs,21.30.-x,21.45.Ff,21.60.De,05.10.Cc 



Introduction. Neutron-rich nuclei are the focus of the ex- 
perimental program of current and next-generation rare iso - 
tope facilities. Emerging phenomena such as halos or neu- 
tron skins make such nuclei ideal laboratories to study nuclear 
interactions in delicately tuned scenarios. These phenomena 
are driven by the evolution of shell structure along isotopic 
or isotonic chains, and motivate the use of ab initio many- 
body calculations to provide a description from first princi- 
ples. Such calculations make it possible to confront modem 
nuclear Hamiltonians from chiral effective field theory (EFT) 
|[l[|2l with a wealth of data beyond few -body systems. 

For s- and D-shell nuclei, the ab initio No-Core Shell Model 
(NCSM) isl 0] provides the necessary capabilities for stud- 
ies of isotopic chains, but for medium-mass nuclei this ap- 
proach is not feasible because of its large computational effort. 
Many-body techniques with more modest computational scal- 
ing, such as the Coupled Cluster (CC) method, can be used 
to probe nuclei in the vicinity of shell closures ||5|-|7|l, but are 
usually not applicable for open-shell nuclei far from shell clo- 
sures. For such nuclei, a self-consistent Gor'kov formalism 
was developed recently JH but at the present stage this 
approach is limited to second-order terms in the many-body 
perturbation expansion. 

In this Letter, we describe the extension of the In-Medium 
Similarity Renormalization Group (IM-SRG) framework of 
Refs. m 

lull to open-shell nuclei by means of a multi- 
reference formulation. We use the resulting MR-IM-SRG 
in conjunction with two other many-body approaches, the 
Importance-Truncated No-Core Shell Model (IT-NCSM) and 
the CC method, to perform the first ab initio study of even 
oxygen isotopes with chiral NNH-3N Hamiltonians. 

Formalism. The main tools for the derivation of the MR- 
IM-SRG are the generalized normal-ordering and Wick the- 
orem introduced by Kutzelnigg and Mukherjee |fl2| . Let us 
denote a string of creation and annihilation operators in ten- 
sorial form. 



and expand it in terms of components that are normal-ordered 
with respect to a given reference state |$), which can be a su- 
perposition of Slater determinants in the multi-reference case 
II13I1 . We obtain 
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where we have introduced irreducible one- and two-body den- 
sity matrices A^i) and A(2): 



A2A4 + A3A4 . 



(3) 



The particle rank of the irreducible density matrices is ev- 
ident from the single-particle indices. Generally, up to n- 
body irreducible density matrices A'"-' appear in the expan- 
sion of an n-body operator, which are defined recursively in 
terms of density matrices of lower rank and encode informa- 
tion about n-body correlations in the reference state II12II . For 
an independent-particle state, all of these matrices except A'^' 
vanish identically. 

Products of normal-ordered operators can be expanded by 
means of a generalized Wick theorem (GWT), e.g.. 
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where = Aj — . In addition to simple contractions 

containing A^^^ and ^'^i'' which also occur in the standard Wick 
theorem, we obtain terms involving A'^^^ . . . , A'-"\ Each den- 
sity matrix must have at least one index from each of the op- 
erators in the product — other terms vanish due to the initial 
normal-ordering ([2]), which greatly simplifies the algebra ifTill . 
In the following, we work in natural orbitals, i.e., the eigenba- 
sis of A^i', where 

^l = niSl, e2^-niSl = -il-m)5l, (5) 

and the eigenvalues are the occupation numbers < Ua < 1. 
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We now consider the IM-SRG operator flow equation 



(6) 



which is used to solve the many-body problem by a continu- 
ous unitary transformation that decouples the ground-state of 
the Hamiltonian H{s) from excitations ifiol [Till . Suppress- 
ing the flow parameter s for brevity, we apply the general- 
ized normal-ordering to H and the generator ?/, and evaluate 
the commutator using the GWT to obtain the zero-, one-, and 
two-body MR-IM-SRG flow equations: 
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where E = ($1 H |$), and the one- and two-body parts of H, 
denoted by / and F, contain in-medium contributions from the 
3N interaction because of the normal ordering lllOi ll ill . The 
symbol [77 ■n- f, T] in Eq. ^ indicates an interchange of the 
one- and two-body parts of rj and H. To close the system of 
flow equations (III-®, we truncate three-body operators 111 ill , 
as well as a term containing A^"^^ in the energy flow equa- 
tion (|2|. We will refer to this truncation as MR-IM-SRG(2). 
Note that the two-body flow equations ^ are identical to the 
closed-shell case ifloHllll . 

As our default choice for the generator, we use the ansatz 
of White I n , Tsl ■ The required matrix elements of the Hamil- 
tonian, such as {^\H : : I''')' which couple the reference 



state to excitations, or ($| '-Af^ ■ H '-All : |<I>), which enter 



the energy denominators, can be evaluated using the general- 
ized normal ordering. This yields 
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where 



G34 = nin2T{2 + nzniVll 

^13 



- {nimT^^ + mniVll + [1 O 2]) . (12) 

The dots in Eqs. JTOl i and ([TTI i indicate terms that are linear 
in A'^'^'. Terms containing higher powers of A(2) or A(") with 
n > 3 are truncated. 

In cases where the flow stalls due to small energy denom- 
inators, we use Wegner's generator 77 = [H, H°'^\ as a fall- 
back, defining the one- and two-body parts of the off-diagonal 
Hamiltonian H°'^ as 
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This generator is free of numerical instabilities but less ef- 
ficient because the flow equations become stiff ifiol [Till . In 
the limit of a single Slater determinant reference state, both 
generators reduce to the forms used for closed-shell nuclei in 

As our reference state, we choose a particle-number pro- 
jected Hartree-Fock-Bogoliubov (HFB) vacuum, \^) — 
PnPz |HFB) lfl6ll . This ansatz for the many -body wave func- 
tion for an open-shell nucleus allows us to enforce spherical 
symmetry in calculations for even nuclei lll7ll . The natural- 
orbital basis of |(1>) is the usual canonical basis of the HFB 
vacuum, allo wing us to use analytic expressions for the den- 
sity matrices llSll . 

We conclude the discussion of the MR-IM-SRG method 
by mentioning that systematic improvements of the truncation 
scheme are possible. One would include 3, . . . , A-body oper- 
ators when expanding the flow equation (|6]l in normal-ordered 
components, as well as additional terms involving irreducible 
density matrices. While the number of flow equations is the 
same as in the single-reference case, their complexity grows 
much more rapidly due to additional terms from the general- 
ized normal ordering 1 13- 12 j . 

Calculation Details. Reference states for the MR-IM- 
SRG calculation are obtained by solving the HFB equations 
in 15 major harmonic-oscillator (HO) shells, and projecting 
the resulting HFB vacuum on good proton and neutron num- 
bers il li I19II . For the 3N interaction, the sum of the HO 
energy quantum numbers of a 3N basis state is limited by 
62 + 63 < i?3,nax = 14, as discusscd in 111 . 
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201. Reduc- 



ing £'3max from 14 to 12 changes the MR-IM-SRG(2) ground- 
state energies for oxygen isotopes by less than 1% for the 
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1. (Color online) Convergence of the MR-IM-SRG(2) ground- 



state energies of O and O with respect to the single-particle basis 
size Cmax, for the NN+3N-full Hamiltonian at Asrg = 2.0 fm~^. 



Hamiltonians used in this work. The intrinsic NNH-3N Hamil- 
tonian is normal-ordered with respect to the reference state, 
and the residual normal-ordered 3N interaction term is dis- 
carded, leading to the normal-ordered two-body approxima- 
tion (N02B), which is found to overestimate oxygen binding 
energies by about 1% 111 , 20ll . 

In this Letter, we use the same nuclear Hamiltonians as in 
our recent IM-SRG and CC studies lEl H^, El] : The NN in- 
teraction is the chiral N'^LO interaction by Entem and Mach- 
leidt, with cutoff Ann = 500 MeV/c il Q. Our standard 
three-body Hamiltonian is a local N^LO 3N interaction with 
initial cutoff A3N = 400 MeV/c. The resolution scale of the 
Hamiltonian is lowered to Asrg = 1.88, ... , 2.2^fm~^ by 



means of an SRG evolution in three-body space 112311 . Hamil- 
tonians which only contain SRG-induced 3N forces are re- 
ferred to as NNH-3N-induced, those also containing an initial 
3N interaction as NNH-3N-full. 

In Fig. [U we illustrate the convergence of the MR-IM- 
SRG(2) ground-state energies for and ^^O with respect to 
the single-particle basis size. At the optimal hil, the change 
in the ground-state energy is of the order of 0. 1 % when we in- 
crease the basis from e,nax — 12 to 14. This rapid convergence 
is representative for all Hamiltonians used in this work. 

Results. In Fig. |2l we show MR-IM-SRG(2) ground-state 
energies of the even oxygen isotopes for NNH-3N-full Hamil- 
tonians with initial cutoffs A3N = 350, 400 and 450 MeV/c. 
For the 3N low-energy constants, we use a fixed cj^ = —0.2, 
and ce = 0.205, 0.098, and —0.016, respectively, which are 
fit to the ''He binding energy in NCSM calculations lEllEell . 
For the NNH-3N-full Hamiltonian with A3N = 400 MeV/c, 
we achieve an excellent reproduction of experimental data all 
the way to the neutron drip line at ^^O 12511, with deviations 
of 1-2%. A recent experiment places the ^^O ground-state 
resonance at < 150 keV above the ^"'O ground-state en- 
ergy I27II . We slightly overestimate this energy in our calcu- 
lation because the HO basis expansion of our single-particle 
wave functions is ill-suited to the description of resonances 
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FIG. 2. (Color online) Dependence of the MR-IM-SRG(2) oxygen 
ground-state energies for the NN+3N-full Hamiltonian on the res- 
olution scale and the initial cutoff Asn. For each A3N, the band is 
obtained by varying Asrg from 2.24 (open symbols) to 1.88 fm~^ 
(closed symbols). Experimental values are indicated by black bars 



and other continuum states. The inset in Fig. |2] illustrates that 
the correct drip-line systematics is independent of Asrg in the 
studied range. The drip line is also robust against variations of 
the cutoff A3N. This suggests that the long-range part of the 
two-pion exchange (2PE) 3N interaction, which remains un- 
changed as we lower A3N, is key to obtaining the proper iso- 
topic trends. The 2PE contribution has significant spin-orbit 
and tensor terms, and is therefore important for the evolution 
of the shell structure along the isotopic chain, as also demon- 
strated in other studies, e.g. 112 811 . 

Let us now discuss the effect of varying the resolution scale. 
,|2C 



As discussed in 111 lLl20ll . the Asrg -dependence of our energies 
is the net result of omitted induced 4N interactions, the Ssmax 
cut, and the MR-IM-SRG(2) truncation of the many-body ex- 
pansion, while the effect of the N02B approximation is found 
to be independent of Asrg- 

For A3N = 350 MeV/c we do not expect significant in- 
duced 4N interactions ll26ll . As Asrg is reduced, we cap- 
ture additional repulsive 3N strength in matrix elements with 
61 + 62 + 63 < i?3max- We also speed up the convergence 
of the many-body expansion and reduce the error due to the 
]VIR-IM-SRG(2) truncation, but for the resolution scales con- 
sidered here, this effect is already saturated. In total, we find 
a slight artificial increase of the ground-state energies as we 
lower Asrg lEitl- 

For our standard choice A3N = 400 MeV /c, effects from 
omitted 4N interactions, the i?3max cut, and the many-body 
truncation cancel, and the AsRc-dependence of the energies 
in Fig. |2]is extremely weak ll ill . The omission of 4N in- 
teractions becomes the dominant source of uncertainty as 
we increase A3N to 450 MeV / c, resulting in an enhanced 
Asrg -dependence of the ground-state energies of the sd-shell 
oxygen isotopes. This observation is consistent with the 



4 



-80- 



-100 



(a) 



NN+SN-ind. 




(b) 



■e— e—e- 



NN+SN-full 




2 4 



6 8 10 12 14 16 18 

A'raax 



2 4 



6 8 10 12 14 16 18 



FIG. 3. (Color online) Extrapolation of IT-NCSM ground-state en- 
ergies for the even oxygen isotopes, obtained using the NN-I-3N- 
induced (left) and NN-l-3N-full Hamiltonians at Asrg = 1-88 fm"^. 
The solid lines indicate the energy extrapolation based on A'„,ax ~ 
8—12 data, while dotted lines guide the eye for smaller A'„,ax- Un- 
certainties due to the importance truncation are smaller than the sym- 
bols used to represent the data. All energies are obtained at optimal 
HQ. 



even stronger Asrg -dependence for A3N = 500 MeV/c, and 
led us to consider the modified 3N interaction with A3N — 
400 Me V /c in the first place lEll S H . 

To assess the quality of our MR-IM-SRG(2) ground-state 
energies, we compare them to results from the IT-NCSM, 
which yields the exact NCSM results within quantified un- 
certainties from the importance truncation [1301 . In the IT- 
NCSM calculations, we use the full 3N interaction without 
N02B approximation, and the iSsmax cut is naturally compat- 
ible with the IT-NCSM model space truncation [jTll. In Fig. [5] 
we show the convergence of the oxygen ground-state energies 
for the NNH-3N-induced and NNH-3N-full Hamiltonians as a 
function of A^max, along with exponential fits which extrap- 
olate iVinax 00 E9U3I. With the exception of ^^O, all 
isotopes converge well, and the uncertainties of the threshold 
and model spaces truncations of the IT-NCSM results are typ- 
ically of the order of 1 MeV. For ^^O, the rate of convergence 
is significantly worse, which is expected due to the resonance 
nature of this ground state. 

In Fig. H we compare the MR-IM-SRG(2) and IT- 
NCSM ground-state energies of the oxygen isotopes, for the 
NNH-3N-induced and NNH-3N-full Hamiltonians with Asrg = 
1.88 fm~^ to experiment. For the latter, the overall agree- 
ment between the two very different many-body approaches 
and experiment is striking: Except for slightly larger devia- 
tions in the drip line nuclei ^^O and ^^O, we reproduce exper- 
imental binding energies within 2-3 MeV. This is a remark- 
able demonstration of the predictive power of current chiral 
NNH-3N Hamiltonians, at least for ground-state energies. For 
further confirmation, we perform CC calculations with sin- 
gles and doubles (CCSD), as well as perturbative triples (A- 
CCSD(T)) miiallllli] for oxygen isotopes with sub-shell 
closures. Using the same Hamiltonians in N02B approxima- 
tion, the MR-IM-SRG energies are bracketed by the CC re- 
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FIG. 4. (Color online) Ground-state energies of the oxygen iso- 
topes from the indicated many-body methods, for the NN-I-3N- 
induced (top) and NN-l-3N-full (bottom) Hamiltonian with A3N = 
400 MeV/c. MR-IM-SRG(2), CCSD, and A-CCSD(T) results are 
obtained at optimal hQ, using 15 major oscillator shells and -Esmax = 
14. The IT-NCSM energies are extrapolations to infinite model 
space. Experimental values are indicated by black bars ll24l[25in . 



suits, and more similar to the A-CCSD(T) values, consistent 
with the closed-shell results discussed in il ill . 

For the NNH-3N-induced calculation, which should be com- 
pared to calculations with the bare chiral NN interaction 
the reproduction of experimental trends fails, and the neutron 
drip line is predicted at wrong mass number, because ^^O is 
bound with respect to ^^O. This illustrates the crucial impor- 
tance of the chiral 3N interaction for a proper description of 
the structure of neutron-rich nuclei ll28ll . 

Let us conclude the discussion by addressing the uncertain- 
ties of our results. The MR-IM-SRG(2) energies lie about 
1.5-2% below the IT-NCSM results. About 1% of this de- 
viation is caused by the N02B approximation. The uncer- 
tainty due to the E^max cut should be smaller than 1 % at low 
Asrg- While the N02B and E'smax uncertainties exhaust the 
greater part of the 1 .5-2% deviation between MR-IM-SRG(2) 
and IT-NCSM, and suggest a very small uncertainty due to the 
many-body truncation, we assume a more conservative many- 
body truncation error of 1-1.5%, and an overall uncertainty 
of our oxygen energies at the level of 3-3.5%, consistent with 
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our closed-shell IM-SRG calculations fllj. Because all ir- 
reducible many-body density matrices vanish in closed-shell 
nuclei, our findings indicate that the truncation of terms con- 
taining A^'^) and higher, or non-linear powers of A'^'^', is negli- 
gible compared to the truncation of induced three-body opera- 
tors. A more detailed analysis of the MR-IM-SRG truncation 
scheme will be presented in a future publication. 

Conclusions. We have generalized the IM-SRG approach 
to multi-reference states, and used the resulting MR-IM-SRG 
method to perform the first ab initio study of the even oxy- 
gen isotopes with chiral NNH-3N Hamiltonians, along with 
the IT-NCSM and the CC method. We achieve remarkable 
agreement with experimental binding energies within theoret- 
ical uncertainties of 3%. This agreement is achieved without 
any re-adjustment of the interaction to experimental data be- 
yond *He, and therefore constitutes an impressive demonstra- 
tion of the predictive power of chiral NNH-3N Hamiltonians. 
The present work also highlights the importance of the 3N in- 
teraction for the nuclear structure of neutron-rich nuclei, as 
demonstrated by the robust reproduction of the oxygen drip 
line. 

While the oxygen isotopes studied here are also within the 
reach of the particle-attached and -removed CC calculations, 
the MR-IM-SRG provides ground-state energies at a signifi- 
cantly reduced computational cost. Thus, it is ideally suited 
for calculations in medium- and eventually heavy-mass open- 
shell nuclei far from shell closures. In such applications, the 
MR-IM-SRG also offers a more complete treatment of many- 
body correlations than other current open-shell methods. 
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